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A Comparative Study of the Parker Instability under Three Models of the 

Galactic Gravity^ 

Jongsoo Kim^'^'^, and S. S. Hong^''^ 
ABSTRACT 

To examine how non-uniform nature of the Galactic gravity might affect length 
' and time scales of the Parker instability, we took three models of gravity, including 

0^ ■ the usual uniform one. In a linear model we let the acceleration perpendicular to the 

' Galactic plane increase linearly with vertical distance z from the mid-plane. As a more 

^ ^ . realistic choice, we let a hyperbolic tangent function of z describe the observationally 

known variation of the vertical acceleration. To make comparisons of the three gravity 
models on a common basis, we first fixed the ratio of magnetic pressure to gas pressure 
I at a = 0.25, that of cosmic-ray pressure at /? = 0.4, and the rms velocity of interstellar 

' clouds at Os = 6.4 km s~^, and then adjusted parameters of the gravity models in such 

' a way that the resulting density scale heights for the three models may all have the 

■ same value of 160 pc. 

. In the initial equilibrium state, the vertical density structure is given by an 

. exponential, Gaussian, and power of hyperbolic cosine functions of z for the uniform, 

i linear, and realistic gravity models, respectively. Performing linear stability analyses 

Qh| onto these equilibria with the same ISM conditions specified by the above a, (3, and 

O ' values, we calculate the maximum growth rate and corresponding length scale for each 

' of the gravity models. Under the uniform gravity the Parker instability has the growth 

■ time of 1.2x10^ years and the length scale of 1.6 kpc for symmetric mode. Under 
^ ■ the realistic gravity it grows in 1.8x10^ years for both symmetric and antisymmetric 

^ . modes, and develops density condensations at intervals of 400 pc for the symmetric 

mode and 200 pc for the antisymmetric one. A simple change of the gravity model has 
thus reduced the growth time by almost an order of magnitude and its length scale by 
factors of four to eight. These results suggest that an onset of the Parker instability in 
the ISM may not necessarily be confined to the regions of high a and /3. 

Subject headings: instabilities — ISM: clouds — ISM: magnetic fields — 
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1. INTRODUCTION 

By performing linear stability analysis upon an initial equilibrium state of interstellar matter 
(ISM) that is supported against an externally given uniform gravity by the pressures of interstellar 
gas, magnetic fields, and cosmic ray particles, Parker(1966) proved that such a system is unstable 
to long wavelength perturbations along the direction of initial unperturbed magnetic fields. 
Since then many studies investigated effects on the Parker instability of rotation (Shu 1974; 
Zwcibel &: Kulsrud 1975; Foglizzo &: Tagger 1994), magnetic microturbulcncc (Zweibel & Kulsrud 
1975), skewed configuration of magnetic fields (Hanawa et al. 1992a), non- uniform nature of the 
externally given gravity (Horiuchi et al. 1988; Giz &; Shu 1993; Kim et al. 1997), self-gravity of 
the ISM (Elmegreen 1982; Nakamura et al. 1991; Hanawa et al. 1992b), and the Galactic corona 
(Kamaya et al. 1997). The Parker instability is expected to play significant roles in the dynamo 
action of accretion disks (Tout & Pringle 1992) and the ejection of mass to galactic halos (Kamaya 
et al. 1996). The instability has been thought an important formation mechanism of the giant 
molecular cloud complexes (GMCs) in the Galaxy (Appenzeller 1974; Mouschovias et al. 1974; 
Blitz & Shu 1980; Shibata & Matsumoto 1991; Handa et al. 1992; Gomez de Castro & Pudritz 
1992). 

As a formation mechanism of the GMCs, the Parker instability under the uniform gravity 
is facing two severe problems. The fastest growing mode has an infinite wavenumber along the 
radial direction, which is the direction perpendicular to both the unperturbed magnetic field and 
the externally given gravity (Parker 1967). Therefore, small-scale chaotic structures rather than 
large-scale condensations are likely to form through the instability (Assco et al. 1980; Kim et 
al. 1998). Another problem lies in the time and length scales. If perturbations of finite vertical 
wavelength are given, at canonical ISM conditions (c/. Spitzer 1978), the Parker instability under 
the uniform gravity grows in a time scale of 1.2x10^ years with the corresponding length scale 
being 1.6 kpc. These scales are too long and too large for the Parker instability to be the formation 
mechanism of the GMCs, since lifetime of interstellar clouds is only about 3x10'^ years (Blitz &; 
Shu 1980) and mean separation of the GMCs is observed to be about 0.5 kpc (Blitz 1991). 

In relation to the second problem of scales, Mouschovias et al. (1974) pointed out that 
physical conditions behind the Galactic shocks are more favorable to trigger the Parker instability 
than in general interstellar space. It is true that in the shocked region one may reduce both 
scales by significant factors. But at the same time an increase in density there makes heretofore 
ignored self-gravity important. The Jeans instability would then override the Parker instability 
(Elmegreen 1982). 

The first problem of an infinite wavenumber has its root at the Rayleigh-Taylor instability, 
whose growth rate increases with increasing wavenumber. Hanawa et al. (1992a) sought a solution 
to this problem of infinite wavenumber from a skewed magnetic field. If the field lines change 
their directions systematically with height from the Galactic mid-plane, the buoyancy is likely 
to be suppressed. According to their analysis, the fastest growing mode of the Parker instability 
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takes a wavelength of A c±10i?, with H being the density scale height. As long as the skewness 
is appreciable, for example, more than 30° /H, the wavelength of maximum growth rate doesn't 
seem to depend sensitively on the degree of skewness. Observational confirmations are still needed 
on the skewness of the Galactic magnetic fields. Even if a significant skewness is confirmed, the 
resulting scale length of lOH is again too large for the GMCs. This brings us back to the second 
problem of scales involved in the perturbations along the field line. 

Parker (1966) assumed, in his pioneering study, the vertical gravitational acceleration to be a 
constant of z. In many ensuing studies the assumption of uniform gravity was almost exclusively 
employed in estimating the time and length scales of the instability. The values of 1.2x10^ years 
and 1.6 kpc we quoted above and A c^lOH of Hanawa et al. (1992a) are all based on the same 
assumption of uniform gravity. However, we do know that the Galactic gravity is of non-uniform 
nature (Oort 1965; Bienayme, Robin &; Creze 1987, hereafter BRC). It is then interesting to see 
how much reductions one may achieve in the length and time scales by introducing a realistic 
model for the Galactic gravity. 

The uniform gravity assumption has brought a few unphysical features to the classical 
picture of the Parker instability. Under a constant acceleration, the equilibrium distribution 
of density follows an exponential function of height z and has a cusp at the mid-plane, which 
accompanies a discontinuity in the pressure gradient at z = 0. Because of the discontinuity the 
linear stability analyses based on the uniform gravity ought to be limited to the perturbations 
of mirror symmetric mode only; while the perturbations of antisymmetric mode are expected to 
dominate the symmetric ones (Horiuchi et al. 1988; Basu et al. 1997). The antisymmetric mode 
would not only slightly increase the growth rate over the value of mirror symmetry, but also 
obviously reduce the distance between condensations by a factor of two. Therefore, models other 
than the uniform gravity are required to make better estimations of the time and length scales for 
the Parker instability in the Galaxy. 

In their linear stability analysis, Giz &: Shu(1993) introduced a hyperbolic tangent function 
of z to model the vertical variation of the Galactic gravity. And recently Kim et al. (1997) used 
a linear function of z as a non-uniform model of the Galactic gravity. Both studies demonstrated 
that dynamical evolution of the perturbations given to the Galactic disk would follow, under such 
non-uniform gravities, two different families of solutions. They are called continuum and discrete 
families. The familiar continuum family has two modes of solutions: One is the original Parker 
mode, which is a slow magnetohydrodynamic (MHD) mode modulated by the gravity, and the 
other is a stable Alfvcn mode. The newly discovered discrete family has three modes: stable fast 
MHD, stable slow MHD, and unstable slow MHD modes. Both studies put their emphases on the 
discovery and characterization of the discrete family, and didn't fully address the question of GMC 
formation in the context of non-uniform gravity. When an effective adiabatic index of the ISM is 
larger than a certain critical value, only the discrete family can have unstable solutions (Giz &; 
Shu 1993; Kim et al 1997). However, the growth time of the unstable discrete solution turned out 
to be about 10 times longer than that of the unstable continuum solution. It is, therefore, difficult 
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to explain formation of the GMCs by the discrete family of solutions. 

In the present study we will concentrate on the continuum family of solutions. Since the 
external gravity is the driving force of the Parker instability, its time and length scales should 
depend on the nature of the Galactic gravity. This line of reasonings motivated us to seek a 
solution to the second problem of scales from the non-uniform nature of the Galactic gravity. 

This paper is organized as follows. In §2 we first introduce three models for the Galactic 
gravity, and then construct equilibrium configuration for each of the gravity models. The MHD 
equations arc linearized in the same section. In §3 dispersion relations arc derived for the three 
gravity models, and the resulting sets of time and length scales are compared with each other. In 
§4 we summarize the paper with some discussions. 



FORMULATION 



Complete MHD equations for the magnetized gas and cosmic-ray particles under an externally 
given gravitational field are given by 



dv 
di 



-V{p + P + 



Stt 



+ —B-VB + pg, 



dB 
'dt 



= V X {v X B), 



dp ^ 

— + • Vp + 7pV ■ V = 0, 
at 

B-VP = 0. 



(1) 
(2) 

(3) 

(4) 
(5) 



Here gas and cosmic-ray pressures are denoted by p and P, respectively. Other symbols have 
their usual meanings. Since the dispersion of random thermal velocities is much smaller than that 
of random cloud velocities, the main source of gas pressure is from the macroscopic turbulent 
motion of the clouds. To describe the local behavior of the Parker instability in the Galactic disk 
we introduce Cartesian coordinates x, y, and z, whose axes are taken to be parallel to the radial, 
azimuthal and vertical directions, respectively. The gravity has only vertical component, whose 
magnitude varies with height z, i.e. g = [0,0, —g{z)]. 

If 7 is larger than one and smaller than a critical value, the system becomes unstable by both 
families of continuum and discrete; if 7 is less than one, there exists convective instability (Kim et 
al. 1997). Since our principal concern is to see the effects of gravity on the Parker instability, we 
simply choose 7 = 1 and limit ourselves to the unstable continuum family of solutions. 
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2.1. Three Models for the Galactic Gravity 

In the Galaxy the main source of gravity is stars, and the vertical gravity can be represented 



.<v'^> . / z 



by 

9r{z) = 2^j^tanhl^^J, (6) 

where <v1> and Ht^ are the velocity dispersion and the scale height of the stars, respectively (Kim 
1990; Giz &; Shu 1993). We call this a realistic gravity as Giz &; Shu did. In addition to this, we 
consider a uniform gravity, 

gu{z)=goi^y (7) 

and a linear gravity, 

gi{z)=g'z, (8) 
where both go and g' are positive constants. 

The uniform, linear and realistic gravity models are represented in Figure 1 by the dotted, 
dashed, and long-dashed lines, respectively. The two solid lines in the figure represent the strength 
of the gravitational accelerations of Oort (1965) and of BRC at solar neighborhood as functions of 
z. The gravity inferred from the distribution of K giants (Oort 1965) increases almost linearly up 
to z ~ 500 pc, beyond which it stays more-or-less constant. The same trend can be found from 
the gravity based on a galaxy model (BRC). 



2.2. Unperturbed States 



Let us suppose that initially an infinitely extended disk of gas and cosmic-ray particles is 
under the influences of magnetic and gravitational fields. The unperturbed magnetic field Bo has 
only an azimuthal component, whose magnitude varies with z, i.e. Bo = [0, Bo{z),0]. Then the 
magnetohydrostatic equilibrium of the system is governed by 

B'oizy 



dz 



Po{z)+Poiz) + 



Stt 



= -po{z)g{z). 



(9) 



To close the differential equation, we take an isothermal equation of state for the gas pressure 
Po = a^Poiz), where the velocity dispersion of clouds is assumed constant everywhere. We also 
assume that the ratio a = Bo/Sirpo of magnetic pressure to gas pressure and the ratio /? = Po/Po 
of cosmic-ray pressure to gas pressure are constants. Then, we have the same z-dependence for the 
unperturbed distributions of density, gas pressure, cosmic-ray pressure, and magnetic pressure as 

Po{z) _ Po{z) _ Po{z) _ Bliz) r Mz) 



where ^{z) is defined by 



PoiQ) Po(0) 52(0) 



exp 



(l + a + /3)a2 



^{z) = [ g{z)dz, 
Jo 



(10) 



(11) 
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and Po(0), Po(0), -Po(O), and Bo{0) denote their mid-plane values. An effective scale height H of 
the density distribution is defined by 



2H 



1 



Po{0) 



Po{z)dz. 



(12) 



Under the uniform gravity, the unperturbed state is described by an exponential function, 



Po,u{^) _ Po,u[z) _ Po,u{z) 



BUz) 



Po,«(0) PoAO) Po,uiO) ^o,«(0) 

with the scale height being 

H(=Hu 



exp 



\z\ 



(1 + a + P)a 



9o 

Under the linear gravity, the unperturbed state is given by a Gaussian function. 



Po,m Po,m Bi^io) 

where the scale height now becomes 

"7r(l + a + /3) 



exp 



'I 



(13) 



(14) 



(15) 



29' 



Hi) = 

For the case of the realistic gravity, it is described by 

Po,r{z) _ Po,r{z) _ fo,r(^) _ -^o,r(^) 
PoAO) ~ Po,r(0) ~ Po,r(0) ~ B^iO) 



sech^^ I — 



z 



(16) 



(17) 



where s is defined by 



{l + a + l3)al' 

If we call {1 + a + I3)a1 as an effective velocity dispersion of clouds, then s is the ratio of the 
velocity dispersion of stars to the effective velocity dispersion of clouds. Substituting equation (IT 



(18) 



into equation (12), we may describe the scale height of clouds in terms of ff* and s, 

(2s 



(2s 



Dir*- h ' 



(19) 



where h is the ratio of the scale height of stars to that of clouds. 

In order to specify the gravity models one should fix all the parameters go, g', H^, and <f^>. 
Under each model of gravity, the magnetized gas and cosmic-ray particles adjust themselves to the 
equilibrium state. This means that the gravity parameters are related to the ISM parameters, a, 
f3, H, and (see eq. [14] for uniform gravity, eq. [16] for linear gravity, and eqs. [18] and [19] for 
realistic gravity). In practice, however, it is difficult to uniquely determine the values of Qq and g' 
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from either Oort's (1965) or BRC's gravity. Furthermore, the values of and <vl> for one type 
of stars are different from those for other type of stars (Mihalas & Binney 1981). Therefore we 
decide to pin down the parameter values on the basis of the ISM conditions. They are specified in 
such a way that the gas scale height resulted from each of the three gravity models may all take 
the same value 160 pc that is known from observations (Falgarone &; Lequeux 1973). In this way 
the comparison of the length and time scales among the three gravity models can be done on a 
common ground. 

We take a = 0.25 and /? = 0.4, which are canonical values of the ISM (Spitzer 1978). For the 
scale height and the rms velocity of interstellar clouds we take 160 pc and 6.4 km s~^, respectively 
(Falgarone & Lequeux 1973). Then, go and g' are equal to 1.4x10^^ cm s~^, and 13.4x10"^ cm 
s~^ kpc~^. In addition to the above values of a, (3, H, and a^, we should also fix if* and <f^> 
according to the chosen value of parameter s. This completes the specification of the realistic 
gravity model. The resulting three models are compared in Figure 1. 

The parameter s defined by equation (18) is exactly the same as R in Giz & Shu (1993). They 
took R fs 3.5 by fixing the value of i?eq ~ 2, ratio of the equivalent half-thickness of stellar disk to 
gas disk. But the equivalent half-thicknesses for various stellar objects in the Galaxy are different 
from type to type. In this paper, however, we will fix the s value in such a way that the resulting 
model of gravity may closely resemble the gravities of Oort (1965) and BRC. As can be seen from 
Figure 1, their gravities are reproduced by the realistic model with s = 16 and s = 9, respectively. 



2.3. Linearized Perturbation Equations 

We limit ourselves to the perturbations that are propagating in the plane defined by the 
azimuthal and vertical directions. In the two dimensional geometry the magnetic vector potential 
A = A(y, z)x is more convenient to use than the magnetic field B, because one scalar quantity, 
A{y,z), is enough to specify the y and z components of the field. Taking the advantage we 
combine the linearized equations (1) through (5) into one for 5A: 



dz^ 



6A 




(2a + 7)Q2 + ^2^2 



7(l + a + /3--j^lnp„. 



1 a 



2 7 \dz 

2 ^ d 



7 



dz 



Inpo 



dy^ 



(20) 



M = 0, 



where is an acoustic wave operator defined by 



5*2 



92 
dy 



2 ■ 



(21) 
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If the following function 



5 A = f{z) exp{itLit — ikyu) 



(22) 



is substituted for the perturbation of vector potential, equation (|20D becomes a second order 
ordinary differential equation 



(2a + -f)kl + ^ ( a + 



7 



dz 



In po 



+agky ^ 2a'yky 



(1 + a + /3)(1 + a + /3 - 7) - -07 



d 

3- In Po 

dz 



1 + 

2 



(2a + -f)uj^ - 2a-ya;ky 
' d^ 



dz^ 



dz^ 
d^ 



In Po 



u;'}f{z)i23) 



+ «7 I ^Inpo 1 >f{z) 



0, 



where f{z) is an amplitude function of the perturbation, u is the angular frequency, and ky is the 
wavenumber along the azimuthal direction. 

As an upper boundary condition (BC), one may set / — at z — -^node- The value of ^^node 
can be either finite or infinite. As a lower BC, one may take either f = or df / dz = at z = 0. 
The former generates solutions having mirror symmetry; while the latter does the ones having 
antisymmetry. Since z-distribution of the unperturbed state has a cusp at the mid-plane under 
the uniform gravity (eq. only one condition, / = 0, is applicable to the lower boundary for 

the case of uniform gravity. 

The effective scale height H and the crossing time over one scale height H/ag are taken as 
normalization units of length and time. Dimensionless variables are then defined by 

n = iLoH/as, C = z/H, Vy = kyH. (24) 

The dimensionless vertical wavenumber, which depends on the lower BCs, should be defined by 

(25) 



2t:{H /2z-^odc) foi^ symmetric modes 
2t:{H / Azjyodc) for antisymmetric modes. 



Since we used the same scale height H and the same rms velocity for all the models, it is not 
necessary to distinguish the dimensionless variables from model to model. 



3. TIME AND LENGTH SCALES 

3.1. Uniform Gravity 

In this paper we will re-derive dispersion relations for the Parker instability under the uniform 
gravity, because we want to point out the problems involved in the time and the length scales, 
and because the results from the uniform model comprise a comparison basis. Under the uniform 
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gravity din po{z) / dz and d'^ In po{z) / dz'^ in equation (23) should be replaced by —1/H and 0, 
respectively (see eq. [O]). In terms of the dimensionless variables, it takes the form 



(2a + 7)1^2 + 2a7z^; 



1 



(1 + a + /?)(! + Q + /3 - 7) - -07 



< = 0, 



(26) 



which is what Parker(1966) gave in his Appendix III. Because the coefficients of d^ fu/dQ"^ and fu 
are constants, and because only symmetric BC at (" = should be applied to the model of the 
uniform gravity, we may set fu oc sin{vzC). Then, an equation for the dispersion relation can be 
written 



1 



+ |2a7(i/2 + 



(1 + a + /3)(1 + Q + /? - 7) - -07 



= 0. 



(27) 



The resulting dispersion relations are shown in Figure 2. The perturbation with a finite 
vertical wavelength grows less rapidly than the one with an infinite wavelength. If we take H 
= 160 pc, Vz = 1-0 corresponds to z^ode — 500 pc. Since the scale heights of interstellar clouds 
and inter-cloud gas are 160 pc and 300 pc (Falgarone and Lequeux 1973), respectively, 500 pc 
is a reasonable choice for the nodal point (e.g., Elmegreen 1982). For the perturbation with the 
infinite vertical wavelength, the minimum growth time of our result is 5.5x10'^ years and Parker's 
(1966) estimate is 3x10'' years. Besides detailed values of a and /3 that went into these estimates, 
the time scale of 3~6xl0^ years is a gross under-estimate; under the Galactic environments it 
seems unrealistic to think of perturbations whose vertical wavelength is much larger than the scale 
height of the cloud distribution itself. In the case of Vz = 1-0, the time and length scales become 
1.2x10^ years and 1.6 kpc, respectively. The growth time, 1.2x10^ years, is much longer than the 
lifetime of GMCs, 3x10'' years (Blitz and Shu 1980), and the length scale, 1.6 kpc, is larger than 
the mean separation of the GMCs, 0.5 kpc (Blitz 1991). 

Figure 3 illustrates the detailed dependences of the the minimum growth time upon the 
parameters a and /3. In most region of the (a, /5) plane the growth time turns out to be longer 
than the cloud lifetime. From the illustration we conclude that under the uniform gravity the 
Parker instability may not play any significant roles in the formation of GMCs. The growth time 
becomes comparable to the lifetime of GMCs only in a limited region of the parameter plane, 
where a and /? take unacceptably high values for the general ISM. This is exactly the reason why 
Mouschovias et al. (1974) invoked the dense region of Galactic shocks as an onset place of the 
Parker instability. 



- 10 - 



3.2. Linear Gravity 

Under the linear gravity, din po^i{z) / dz becomes —{tt/2){z/H'^) (see eq. [jl5|). After 
substituting it and its gradient into equation (23), we can write the resulting equation in the 
following dimensionless form: 

^ + {Ei- Vo,i e)fi = 0, (28) 

where Ei and Vo^i are given by 

TT 04 ^ [(2a + + f (1 + g + m'' + 2«7^,' 



4 (2a + 7)1^2 + 2a7t'2 

1 (1 + a + /?)(! + a + /3- 7)1^2 



7r2 



(30) 



(2a + 7)rj2 + 2a7z^2 

To derive the dispersion relation we will use the same method given in the Appendix of Kim et al. 
(1997). 

Kim et al. (1997) showed that, under the linear gravity, the dispersion relation for the 
continuum family of solutions with a mirror symmetric lower BC are nearly the same as that 
with an antisymmetric BC. Under the point-mass-dominated gravity, however, the antisymmetric 
modes grow faster than the symmetric ones (Horiuchi et al. 1988). Such difference in the behavior 
of growth rate stems from the difference in the nature of the chosen gravity models. The growth 
rate under the point-mass-dominated gravity is sensitive to the lower BCs, because the gravity has 
its maximum close to the lower boundary. On the contrary the growth rates under the linear and 
realistic gravities are insensitive to the lower BCs, because these models have their maxima at far 
from the mid-plane. So in this paper we present the dispersion relations only for the symmetric 
lower BC. 

In the linear model the vertical acceleration increases with z without a bound, while in the 
Galaxy it approaches a finite value at large distance from the mid-plane. We should, therefore, 
limit the disk extent within a finite height from the central plane. The dispersion relations shown 
in Figure 4 are for the perturbations with = 0.5 and = 1.0. This figure clearly indicates a 
strong dependence of the growth rate on the vertical wavenumber. We think = 1.0 a realistic 
choice for the Galaxy. With the same parameters that are used in Figure 2, the perturbation of 
maximum growth rate has a horizontal wavelength 340 pc and a minimum growth time 1.5x10'^ 
years. 

The linear gravity drives the Parker instability much faster than the uniform gravity. By 
almost an order of magnitude reduction is achieved in the growth time scale, even with the 
perturbation wavelength as short as 340 pc. Such reduction is possible for a wide range of a and 
[3 values. This can be seen from the contours, in Figure 5, of equal minimum growth times traced 
out in the (a,/3) plane. Comparison of the two sets of figures (Figs. 2 and 3 versus Figs. 4 and 5) 
indicates that the growth rate is indeed sensitive to the nature of the externally given gravitational 
fields. 



3.3. Realistic Gravity 




with 



(G3 + G4)/Gi, 
(G3 - G2)/Gi 



(32) 
(33) 



where Gi, G2, G3, and G4 are given by 



Gi = (2q + 7)1^2 + 2a7i/2 



(34) 



G2 = [-(2 + 2/3 - 7)f]2 + 2a7z.2] 

G3 = {-(2a + 7)1^2 + [4(1 + a + /3)(l + a + /?-7)- 2a7]z.2| 
G4 = -0^ - (2a + -1)^1^^ - 2aju^. 



(35) 



(36) 



(37) 



To evaluate the dispersion relations we again use the method explained in the Appendix of Kim 
et al. (1997). 

The resulting dispersion relations of the symmetric mode are shown in Figure 6 for different 
values of s. The larger the s value is, the more unstable the system becomes. This is because the 
gravity with a larger value of s is stronger than that with a smaller one (see Fig. 1). To show the 
effect of the vertical wavenumber on the growth rate, the dispersion relations for the cases of i^z 
= 0.0, 0.5, and 1.0 are compared to each other in Figure 7. The smaller the vertical wavenumber 
or the longer the vertical wavelength is, the faster the instability grows. For the set of parameters 
s = 9, = l,a = 0.25,/? = 0.4, and 7 = 1, the minimum growth time becomes 1.8x10'' years 
and the corresponding wavelength is 400 pc. If s = 16 is taken with the other parameters being 
fixed at the same values, the time and length scales are slightly reduced. Figure 8 shows the 
equi-growth time contours in the (a, /3) plane for the case of 7 = 1, i/^ = 1, and s = 9. Under the 
realistic gravity the growth time becomes shorter than the GMC lifetime for most range of the a 
and P values. 



The dispersion relations resulting from the three models of gravity are compared in Figure 9 
The same set of parameter values, a = 0.25, /? = 0.4, 7 = 1.0, I'z = 1-0, is used for all the three 
models. To do a good justice in the comparison, the time and length scales of the three models 



3.4. Comparison of Time and Length Scales 
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are normalized by using the same set of effective scale height 160 pc and rnis velocity 6.4 km s . 
The vertical wavenumber Vz = 1 places a nodal point at z ~ 500 pc. Among the three the linear 
gravity gives the strongest acceleration near the nodal point; while the uniform gravity does the 
weakest (see Fig. 1). Since it is the external gravity that drives the Parker instability, the growth 
rates for the linear gravity are generally higher than for the other two. 

Table 1 lists the minimum growth times and their length scales for each of the three gravity 
models. We include the scales for both the symmetric and antisymmetric modes. The parameters 
used in the calculation of these scales are the same as those in Figure 9. The length scale 
corresponds to the inter-distance of the condensations that would be formed by the Parker 
instability. So the length scale for the symmetric mode, which is equal to the perturbation 
wavelength, is twice the value for the antisymmetric one. We should also clarify the vertical 
wavenumber used for the antisymmetric modes. As is done for the symmetric modes, the upper 
nodal point is set at z ~ 500 pc for the antisymmetric ones. Since the antisymmetric mode doesn't 
have a nodal point at z = 0, this choice corresponds to the vertical wavenumber = 0.5, which is 
smaller by a factor of two than that for the symmetric one. It is, however, fair to use the scales 
with the same upper node rather than those with the same vertical wavenumber. This point 
should be kept in mind, when one compares the scales of the time and length for the symmetric 
modes with those for the antisymmetric ones. 

In spite of the incompatibility of the antisymmetric mode with the uniform gravity, we 
have listed, in Table 1, the time and length scales for the antisymmetric mode, because the 
incompatibility stems from an over-simplification of the Galactic gravity, and because we want to 
compare the scales from different models of gravity. Since the realistic gravity model with s = 9 
represents the results of BRC fairly well, we take 1.8x10'' years for the growth time scale and 400 
pc (200 pc) for its length scale of the symmetric (antisymmetric) mode of the Parker instability. 
These time and length scales are about 1/7 and 1/4 of those for the uniform gravity, respectively. 

4. SUMMARY AND DISCUSSION 

To investigate the effect of external gravity on the continuum family solutions of the Parker 
instability, we introduce uniform, linear, and realistic models for the Galactic gravity. The gravity 
models are specified by fixing the ratio a of magnetic pressure to gas pressure at 0.25, the ratio 
P of cosmic-ray pressure to gas pressure at 0.4, the rms velocity < >^/^ of interstellar clouds 
at 6.4 km s~^, and the effective scale height H of cloud distribution at 160 pc. For the realistic 
model based on the hyperbolic tangent function, we need to specify one more parameter s, which 
is the ratio of the velocity dispersion of stars to the effective velocity dispersion of clouds (see eq. 
|l8|| ). The Galactic gravities of Oort (1965) and BRC at solar neighborhood are reproduced by the 
realistic model with s = 16 and s = 9, respectively. 

Under the uniform gravity the time and length scales of symmetric mode solution are 1.2x10*^ 
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years and 1.6 kpc; while they become 1.8x10^ years and 400 pc under the realistic gravity with 

s = 9. Under the linear and realistic gravities the antisymmetric modes grow at almost the same 
rate as the symmetric mode. However, as for the separation of condensations, the length scale of 
the antisymmetric mode becomes, of course, a half of the symmetric mode's. Changing the nature 
of gravity from uniform to realistic has thus reduced the time and length scales in the Galactic 
environments by factors of 7 and 4, respectively. Therefore, it is not necessary to invoke high 
values of a and (3 to reduce the scales. 

The rotation of the Galactic disk is known to exercise a stabilizing effect on the Parker 
instability (Foglizzo &; Tagger 1994). With the uniform gravity Zweibel & Kulsrud (1975) found 
that the rigid body rotation would increase, over the non-rotating case, the growth time by a 
factor of two. With a gravity model based on a combination of linear and hyperbolic tangent 
functions, which is very close to ours, Hanawa et al. (1992b) found that the rigid body rotation 
would increase the growth time by factors of 1.3 to 7.5 for a varying from 1.0 to 0.1. They placed 
the upper boundary of the Galactic disk at z = 200 pc, while we did at 500 pc. Because the 
external gravity drives the gas to slide down the field lines from top part of the disk first, the 
growth rate one obtains from linear stability analysis depends rather sensitively on the magnitude 
of acceleration at the upper boundary. The acceleration at z = 200 pc is about two thirds of its 
value at z = 500 pc. Therefore, an increase in the growth time the rotation would bring to our 
case may not be so large as their findings. Furthermore, they didn't include cosmic-ray particles as 
an ISM constituent. Since the "light" cosmic-rays drive the field lines to buckle mainly upward, a 
triggering of the instability by the cosmic-ray pressure doesn't seem to be seriously affected by the 
Galactic rotation. If the Galactic rotation had been included in our analysis with the non-uniform 
gravities, the growth time would have increased by a factor of 2 or 3 for moderate values of a and 
f3. However, for very low values of a and /3, the rotation is very likely to make the growth time 
longer than the GMC lifetime (Figure 8). 

Are the GMCs formed by the Parker instability? As far as the time and length scales are 
concerned, our linear stability analysis done with the realistic gravity model gives a positive 
answer to the question. To have a better answer we should also know how much enhancement 
in density can be made by the Parker instability. From the linear analysis alone one may not 
have information on the density enhancement factor. Very recently Basu et al. (1996; 1997) 
performed two-dimensional simulations of the Parker instability under the uniform gravity. The 
density enhancement they obtained at the mid-plane amounts to a factor of only 2. Under 
the non- uniform gravity the enhancement factor may not be much larger than 2, because the 
acceleration of any realistic gravity models anyhow approaches zero near the mid-plane. As for 
the density enhancement the answer is likely to be a negative one. The GMCs may have formed 
through a cooperative interplay of the Parker and Jeans instabilities. 
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Table 1. Comparison of Time and Length Scales 



scales 




uniform 




linear 


realistic (s = 9) 




symmetric 


antisymmetric* 


symmetric 


antisymmetric 


symmetric 


antisymmetric 


time [year] 


1.2 X 10® 


1.2 X 10® 


1.5 X 10^ 


1.5 X 10^ 


1.8 X 10^ 


1.8 X 10^ 


length [pc] 


1600 


800 


340 


170 


400 


200 



'Incompatible with the exponential distribution of the unperturbed state 
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Fig. 1. — Vertical acceleration at solar neighborhood. The two solid lines represent the results 
of Oort (1965) and of Bienayme, Robin &: Creze (1987). A dotted line, a dashed line, and three 
long-dashed lines are for the models of uniform, linear, and realistic gravities, respectively. The 
parameter s is the ratio of the velocity dispersion of stars to the effective velocity dispersion of 
clouds. 

Fig. 2. — Dispersion relations of the Parker instability under the uniform gravity. Each curve is 
marked by the value of vertical wavenumber. ly^ = 0.0 corresponds to an infinite wavelength. If 
the density scale height of clouds is 160 pc, = 0.5 and 1/^=1.0 have nodal points at ~ 1000 pc 
and ~ 500 pc, respectively. The ordinate represents square of the normalized growth rate, and the 
abscissa does the normalized horizontal wavenumber. The effective adiabatic index 7, the ratio a 
of magnetic to gas pressure, and the ratio /3 of cosmic-ray to gas pressure are specified within the 
frame. 

Fig. 3. — Loci of equi-growth time are traced out in the {a, (3) plane for the Parker instability under 
the uniform gravity. 

Fig. 4. — Dispersion relations of the Parker instability under the linear gravity. Each curve is 
marked by the vertical wavenumber Uz- The ordinate represents square of the normalized growth 
rate, and the abscissa the normalized horizontal wavenumber. The system parameters used in the 
calculation are given in the frame. 

Fig. 5. — Loci of equi-growth time are traced out in the (a, (5) plane for the Parker instability under 
the linear gravity. 

Fig. 6. — Dispersion relations of the Parker instability under the realistic gravity. Each curve is 
marked by the value of s, which is the ratio of the velocity dispersion of stars to the effective 

velocity dispersion of clouds. The ordinate represents square of the normalized growth rate, and 
the abscissa the normalized horizontal wavenumber. The system parameters are specified within 
the frame. 

Fig. 7. — Dispersion relations of the Parker instability under the realistic gravity. Each curve is 
marked by the vertical wavenumber Vz- The ordinate represents square of the normalized growth 
rate, and the abscissa the normalized horizontal wavenumber. The system parameters are given in 
the frame. 

Fig. 8. — Loci of equi-growth time are traced out in the (a, 13) plane for the Parker instability under 
the realistic gravity. 

Fig. 9. — Dispersion relations of the uniform, linear, and realistic gravities are compared with 
each other. The ordinate represents square of the normalized growth rate, and the abscissa the 
normalized horizontal wavenumber. The system parameters are specified in the frame. 
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